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We propose a probabilistic quantum algorithm that decides whether a monochrome picture 
matches a given template (or one out of a set of templates). As a major advantage to classical 
pattern recognition, the algorithm just requires a few incident photons and is thus suitable for very 
sensitive pictures (similar to the Elitzur-Vaidman problem). Furthermore, for a 2" x 2™ image, 
0(n + m) qubits are sufficient. Using the quantum Fourier transform, it is possible to improve 
the fault tolerance of the quantum algorithm by filtering out small-scale noise in the picture. For 
example images with 512 x 512 pixels, we have numerically simulated the unitary operations in order 
to demonstrate the applicability of the algorithm and to analyze its fault tolerance. 

PACS numbers: 03.67.Lx, 42.50.-p, 42.30.Sy. 



I. INTRODUCTION 

It is well known that quantum computers are suited 
to solving certain classes of problems much better than 
classical computers. A prominent example is Shor's al- 
gorithm for factoring an integer number with an ef- 
fort that grows polynomially in the number of its digits, 
which is believed to be classically impossible. A furthcr 
impressive example is Grovers algorithm Q for finding 
an item in an unsorted database: There, the effort grows 
only as the square-root y/N of the number N of database 
entries on a quantum computer, whereas it grows with 
A" on a classical computer. In addition, there are fur- 
ther black box problems such as the Deutsch kiand the 
Deutsch-Josza algorithm Q as well as others 
an overview see e.g., Q). 

The possible speedup of quantum algorithms is es- 
sentially enabled by the feature of quantum parallelism. 
This parallelism permits to calculate with a superposition 
of states on a quantum computer, which is not possible 
on classical computers. The first quantum computers 
have already been constructed. For example, Shor's al- 
gorithm has been used on an NMR quantum computer 
[9j to factorize the number 15. This is certainly not im- 
pressive if one considers the smallness of the number but 
nevertheless serves as a proof of principle. 

However, there exists a plethora of further classically 
challenging problems such as e.g., pattern recognition 
|lf)j | , which can also benefit from the application of quan- 
tum algorithms 0, 0, ^| . It has already been demon- 
strated that quantum parallelism can be exploited to 
idcntify and localize a regular simple pattern 14] within 
an otherwise unstructured picture. Here, we will present 
a probabilistic quantum algorithm that is capable of rec- 
ognizing an arbitrary image (even in the presence of noise 
at a moderate level) after shining a few photons on the 
picture. In contrast to previous pattern matching ap- 
proaches (see, e.g., 12]), it does not require the copying 
of quantum states (which is only possible probabilisti- 
cally for general states) and enables a (probabilistically) 
non-destructive measurement, cf. 



II. PROBLEM DEFINITION 

Following the problem description of 0], ^ us con ~ 
sider a large rectangular N x x N y array of unit celis that 
may either be black (absorptive) or white (reflective) . To 
allow for a binary representation, we will consider cases 
where the number of array celis in every dimension are 
powers of 2, i.e., log 2 N x / y = n x / y with integers n x / y . 
Further-on, we will denote the pixels in the array that 
are reflective (white) as points. 

The problem consists of recognizing whether the pat- 
tern in the array matches a given template (for exam- 
ple, a letter of an alphabet). Note that this slightly dif- 
fers from template matching discussed e.g., in [Ï2l Í13| . 
where a template is to be found that optimally matches 
the given unknown quantum state. With using varying 
templates however, one can evidently establish a relation 
between these problems. Some aspects of the presented 
quantum algorithm will be similar to the discrimination 
problem of known quantum states |Ï5j via suitably cho- 
sen positive operator valued measures (POVM), which 
will be discussed in section IVni 

The classical approach to measure the displayed pat- 
tern on the array would be to shine light (consisting of 
many, say 0{N X x N y }, photons) on the array and to 
measure absorption and transmission accordingly. How- 
ever, in the case we consider here, the array is also as- 
sumed to be very sensitive (imagine, for example, an ex- 
posed but not yet developed film or a pattern of partially 
fluorescent ions in a Paul trap), such that each absorbed 
photon causes a certain amount of damagc. Evidently, 
the classical measurement approach would significantly 
disturb the system. 

One might worry that the momentum of the reflected 
photon (recoil) also disturbs the picture, but this effect 
will be extremely small if the mass of a pixel is sufficiently 
big. Alternatively, one might imagine the reflective pixels 
to be transparent and to place a mirror behind the array. 

Then, a quantum algorithm can cope with such a task: 
The feature of quantum parallelism can be exploited by 
storing the relevant information of the image in a quan- 
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tum superposition state with using just a single photon 
and not (or only very little) destroying the image. The 
task is thus similar to the Elitzur-Vaidman problem 
which allows for testing for the existence of an object 
without any energy-momentum exchangc. 



III. READ-OUT SCHEME 

Generally, the coordinates (x, y) of a pixel in the image 
can be written in their binary representation 



x = X\ ® . . . ® x nüs 



y = yi® ...®y nv 



(1) 



where Xi and yj denote the i-th and j-th bit of x and 
y, respectively, i.e., x = X\2 nx ~ l + a^" 31-2 + . . . + x nx . 
Regarding these coordinates (x, y) as control qubits 
|ar) = |xi)|x2) . . . \xn x ), the photon probing the image can 
be entangled with the coordinates in the following way: 
The photon passes through a series of quantum con- 
trolled refractors R® x . . . RfRn y ■ ■ ■ R\ that effectively 
displace the photon by defined distances Axí or Ay,, if 
the control qubit \xí) or \yi) is in the state |1) and leave it 
unaffected if the control bit is in the state |0). By choos- 
ing the displacement of the refractor Rf as Axí = 2 l Aa;o 
(and likewise for the other refractors Rf), the final dis- 
placement of the photon corresponds to the coordinates 
(x, y) which encode the position of a single image pixel. 



Ix,> lxp> |x 3 > 





Ix, > lx 2 > IXj> 



FIG. 1: [Color Online] Left: Illustration of a series of three 
quantum-controlled refractors required for an 8-pixel one- 
dimensional array. A photon created by the source (box in 
the left) passes the refractor series, interacts with the pixel 
array and takes - if reflected - the reverse path. Depending 
on the value of the control qubits \xï), the refraction index 
of the medium is changed, which effectively produces a dis- 
placement of the photon. Right: Schematic representation 
of the two-dimensional configuration (exemplified for an 8 x 4 
array). The first three refractors generate the displacement 
in a;-direction (for 8 pixels) and the remaining two act in the 
y-direction (for 4 pixels). 

In a laboratory setup, this could for example be re- 
alized by using a varying refractor thickness, see figure 
n left panel. For quantum controlled devices, one can 
generate superposition states of several pixel positions 
in the image as control qubits. Storing these bits in a 
coherent superposition state containing all pixels with 
equal amplitudes, a single photon can be forced to in- 
teract simultaneously and uniformly with the complete 
array. 



One should be aware that this causes a strong entan- 
glement between the photon position and the control 
qubits. Thus, any measurement on the photon would 
affect the refractor control qubits as well and completely 
decohere them. The interaction of the photon with the 
image corresponds to absorption or reflection. In case 
of (perpendicular) reflection, the photon will reverse its 
original path (this is enforced by the entanglement with 
the refractor control qubits) and after passing through 
all refractors again, all information about the way of the 
photon (given that it was reflected) is lost. Thus, the en- 
tanglement with the control qubits is partially removed 
(quantum eraser) and we obtain a coherent superposi- 
tion of the reflecting pixels (points). In the other case 
(absorption), the entanglement cannot be removed, and 
the algorithm fails - i.e., one has to send another photon 
if permitted by the fragility of the quantum array. As a 
formal simplification, one can consider the action of the 
configuration of the setup in figure a quantum black 
box: 



B : \x) <8> \y) <g> |0) -> \x) <g> \y) <g> \f(x,y)) 



(2) 



which encodes the output in the characteristic function 
f(x, y) of the image. The function f(x, y) takes the value 
1 if the pixel x <£> y is reflective (i.e., if \x) ® \y) encodes 
a point) and if the pixel is black (i.e., if the photon 
is absorbed). If the control qubits on the refractors are 
initially prepared in a superposition state (e.g., acting 
Hadamard gates on each qubit), the characteristic func- 
tion f(x,y) of the image is tested for all pixels simulta- 
neously 
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H^lot»»)) ® |0< w »>) ® |0) 

N X -Í Ny-Í 



-T^W £ È \x)®\y)®\f{x,y)). (3) 



Measuring the third register (i.e., the existence of a re- 
flected photon) and obtaining |1) as a result prepares the 
quantum state as an uniform superposition of all points 
in the image. The other outeome (|0)) corresponds to 
the absorption of the photon and would lead to entan- 
glement between the refractor control qubits and the 
image. Thus, with the outeome |1) (outgoing photon) 
one has prepared a quantum state (in the refractor con- 
trol qubits) that is suitable for performing further cal- 
culations. An alternative scheme based on linear òptics, 
which could be applied if the image is not too large, is 
discussed in the Appendix. Note that this scheme has 
the advantage that the image does not have to be loaded 
into a possibly fragile quantum memory (cf. [Ï2l Il3j|). 
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IV. QUANTUM ALGORITHM 

After extracting the superposition state containing all 
points of the image 
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(4) 



where M with 1 <C M < N x N y denotes the total number 
of points, one has to decide whether it corresponds to a 
given template. 

Obviously, a N x x N y array (in black and white) could 
contain 2 NxN " different images - but the Hilbert space 
of all possible quantum state has merely N x N y = 2 ni+rly 
dimensions. Hence different images will not correspond 
to orthogonal quantum states in general and thus it is 
- even in the absence of noise - not possible to distin- 
guish thcm with certainty, i.e., the presented quantum 
algorithm can only work probabilistically. 

The problem of efhciently discriminating non- 
orthogonal quantum states has been studied exten- 
sively and is usually formulated within the framework 
of positive operator valued measures (POVM): Given 
a quantum state out of a set of perfectly known (non- 
orthogonal) quantum states, it is possible to construct a 
POVM which decides probabilistically which one of these 
states it is with minimum error probability, see section 
IVIII However, the design of a suitable POVM will be - 
especially in the case of large images (high-dimensional 
Hilbert-spaces) with more than two possible quantum 
states - quite demanding in general |15|, [l6j. In addi- 
tion, if the quantum states are not known exactly - e.g., 
in the presence of perturbations and noise - the construc- 
tion of the optimal POVM is even more complicated. Fi- 
nally, the physical realization of such an optimal POVM 
will also be quite demanding in general and can easily 
destroy part of the speed-up of the quantum algorithm. 
Therefore, we propose a different procedure that shifts 
the effort from designing a suitable POVM towards the 
implcmentation of a quantum oracle gate (in an inverse 
application of Grovers algorithm). 

Using the characteristic function f(x,y) of the image 
(which yields 1 for the marked items, i.e., points of the 
image, and otherwise) as an oracle, a Grover iteration 
generates a rotation in the two-dimensional sub-space of 
the Hilbert space spanned by the coherent superpositions 
of all numbers \s) = H {n ^ <g> \0^) on the 
one hand and the numbers of the marked items, i.e., the 
state \ in the above equation, on the other hand. Since 
there are M solutions to this search problem, 
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(5) 



inverse Grover iterations will rotate the original state | \E r ) 
in Eq. (@J into the total superposition state \s), or at 
least close to it (discretization error, cf. After ap- 

plication of a Hadamard gate to each of the qubits, the 



resulting state would be lO^™*)) £g> |0( n ")). If the image 
indeed perfectly equals the template - i.e., if thcir char- 
acteristic functions coincide - and if the Grover iterations 
were assumed to be perfect (no discretization error), the 
state | \E f ) in Eq. (JIJ prepared after measuring the return- 
ing photon would be unitarily transformed into the final 
state 

| («*)) ^ |0("v)). Consequently, a measurement in 
the computational basis would yield zeros for all bits. 

If the image differs from the template (and assuming 
a perfect Grover rotation E|), the probability of ending 
up in the final state |0^™ x ') (£> |0' raï · ) ) would be given by 
the overlap of the characteristic functions of the template 
f(x,y) and the image f'(x, y) 
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x=0 y=0 



(6) 



since all involved operations were unitary and hence 
probability conserving. Since this overlap determines 
how similar image and template are, measuring all qubits 
at the end and obtaining zero everywhere is a strong in- 
dication that the image equals the template or at least is 
very similar to it. Obtaining 1 somewhere, on the other 
hand, indicates that the image does not equal the tem- 
plate with high probability |l9j . As usual, the result can 
be made more decisive by repeating the whole algorithm. 



V. NOISE FILTERING 

So far, the input state l'í) was assumed to be perfect. 
However, in reality neither the realization of the image in 
the array nor the reflection and absorption properties of 
the image pixels can be assumed as perfect. Therefore, 
the bàsic algorithm described above can be improved us- 
ing the quantum Fourier transform to reduce possible 
perturbations: As in classical pattern recognition, one 
can perform a cutoff of large wave-numbers, which re- 
duces noise such as pixel defects. In a quantum algo- 
rithm, such a cutoff has to be implemented in a suitable 
way - just measuring the fc-value and subsequently de- 
ciding whether it lies above of below the cutoff does not 
work, since a hard fc-measurement completely decoheres 
the quantum state. Therefore, a noise-filter should be 
realized via the measurement of an ancilla qubit which 
has been coupled to the control qubits: To this end, we 
introduce a unitary noise-filter operator 

Af[\k x )<E> \k y }® |0)] = 
\k x ) ® \k y ) ® {cos[tf(k x ,k y )}\l) + sm[$(k x ,k y )}\0}} , (7) 

with a suitably chosen noise-filter function $(k x , k v ). 
Measuring the last register and obtaining |1) performs 
the desired cutoff in frequeney space. 

Let Mtp denote the number of template points and 
M{ m w M tp the number of image points. Then, the quan- 
tum algorithm can be summarized as follows (compare 
also figure[21): 
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1. initialize the state vector with uniform superposi- 
tions in x and y as well as two ancilla qubits 

I*!) =n (nx) \o n -} ®wK)|o™«) ® |o) ® |o) 

2. apply the black box with the first ancilla qubit 
|* 2 ) =B|*i) 

3. measure the first ancilla qubit and proceed after 
obtaining the outcome |1) 



I*.,) = M; 



im 1/2 E \x) ® |v) ® |i) ® |o) 

(a:,y)£image 



4. perform quantum Fourier transform in |x) and |y) 



|*4> = QFTxJVs) 



5. apply noisc filter operator with second ancilla 
|* 5 ) = A^ 4 ) 

6. measure the second ancilla qubit and proceed after 
obtaining the outcome |1) 

|i> 



Efc T .fc„ /cutC^xj fcj/)|fcx) 



fc?/ 



1 



1 



with f cut (k x ,k y ) 



7. perform the inverse quantum Fourier transform 



|* 7 ) = QFTl y \* 6/ 



i. perform R 



l*8> =a- fi i*7) 



A/t, 



inverse Grover iterations 



9. apply Hadamard gates on non-ancilla qubits 

|* 9 )=W (fl " ) «l« (n,) |*8) 

10. measure final state in computational basis 
|* 9 ) I |0("-)) ® |0("f)) <g> |1) ® |1) 

The algorithm fails when measurement of one of the an- 
cilla qubits yields |0), i.e., if the photon is absorbed or if 
a projection onto the wrong fc-values is performed. 
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FIG. 2: [Color Online] Quantum circuit of the template 
matching algorithm for 14 control qubits (7 for x and y re- 
spectively) plus 2 ancilla qubits (ancilla qubit vàlues shown 
for a successful run). 

The probability of obtaining finding the system in the 
state |0 . . . 0) ® |1) (g> |1) as the result of the last measure- 
ment corresponds to the probability of template match- 
ing, i.e., given a template such as the letter A and 



an initial state such as the (noise-perturbed) state \B), 
the quantum al gori thm decides whether the template is 
matched or not [1911. 



VI. ALGORITHMIC PERFORMANCE 

Starting from a possibly noise-perturbed initial state 
as prepared by step 3 of the quantum algorithm, we 
have numerically simulated the action of the correspond- 
ing unitary gates and measurements for the 512 x 512 
example arrays in figurc [3] As the number of allowcd 
states increases exponentially with the number of sim- 
ulated n x + n y qubits, the numerical simulations would 
involve 2™ x+ ™« x 2 ïla:+rl « matrices which do not fit into 
the main memory. Fortunately, the involved unitary op- 
erations can be expanded into combinations of one or two 
qubit operations 8], which can be calculated. Thus, the 
whole algorithm for pattern recognition for n x + n y = 18 
qubits runs in a time of few seconds. As a noisc- filter 
function a sharp cutoff was used, i.e., 



${k x , ky) — 




n 

othcrwisc 



(8) 



which leads to a simple projection on the allowed fc- 
values. In addition to the high-frequeney (fc > fc ma x) 
components, it can be advantageous to remove the 



component as well, especially if the noise 



significantly changes the total number of points in the 
image. 

The computational complexity of the algorithm de- 
pends on the total number of pixels and on the num- 
ber of points in the template. The Hadamard gates and 
the quantum controllcd refractors require 0{n x + n y } 
operations and the quantum Fourier transforms in- 
volve 0{n x +n y } gates. The necessary number of 
Grover iterations depends on the size of the template 
R = 0{ y/N x N y /M t p}, and the number of involved gates 
per Grover iteration depends on the physical realization 
of the oracle function f(x,y). Many template points 
(e.g., bold and large letters) lower the number of required 
Grover iterations. 

Similarly, one has to estimate the failure probability, 
i.e., how often measurement of one of the ancilla bits 
yields zero. The probability for the photons to be re- 
flected (first ancilla) is given by the ratio of the image 
size (number of points) over the total area of the array 
(number of pixels), i.e., many image points are favor- 
able (similar to the above argument). For the second 
ancilla the failure probability depends on the amplitudes 
of the removed fc-values and thus roughly on the amount 
of noise and the number of fine details in the image. 

For specific templates we have made our algorithm ex- 
plícit: With a given template A or B (cf. figure top 
panels), we numerically calculated the recognition prob- 
ability for all the possible combinations (image A and 
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template B versus template A ctc.) for various noise 
levels (cf. figure lower panels). 



decreases much faster. 
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FIG. 3: [Color Online] Left: Images of the input states 
(512 x 512) with various noise levels. From top to bottom 
the image pixels of the initial state have been inverted with 
a probability of 0%, 5%, 10%, 20%, and 40%, respectively. 
Right: Panels show the distribution of the squared amplitude 
of the inverse of the Fourier-smoothed images (i.e., after step 
6 of the quantum algorithm). Only k- vàlues with < k < 40 
have been kept. All amplitudes larger than half the maximum 
amplitude in every image are shown in red, whereas the thin 
black isoline encodes the quarter of the maximum amplitude. 

The results of the numerical simulations are shown in 
figure01 The slight differences between the letters A and 
B result from the different numbers of image points occu- 
pied by the two templates (discretization error [lflj of the 
Grover iterations etc.), as has been checked by applying 
the algorithm to a different (more symmctric) font. An 
important quantity is the the discrimination ability of 
the algorithm, i.e., the capability to distinguish between 
the alternatives in the given alphabet A and B, see fig- 
ure 21 A rough measure for the discrimination ability 
is given by the difference between the acceptance proba- 
bility when the input state corresponds to the template 
(circle symbols) and the acceptance probability when the 
input state does not match the template (squarc sym- 
bols). Without noise, this difference is over 70% and 
even with 40% noise, the discrimination ability is still 
significant (around 30%), provided that a noise filter is 
applied. Without noise filter, the discrimination ability 
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distortion of given letter 

FIG. 4: [Color Online] Probability to recognize the initial 
state as the template versus noise level. Small hollow symbols 
represent results without noise-filter applied. 



VII. FAILURE PROBABILITY 

As already mentioned in sections |n] and IIVI the task 
in section IIVI is similar to the well-known discrimina- 
tion problem: Given a quantum state out of a set of 
perfectly known (non-orthogonal) quantum states, it is 
possible to construct a positive operator valued measure 
(POVM) which decides which one of these states it is. 
Since non-orthogonal states cannot be distinguishcd with 
ccrtainty, such a POVM measurement can never be per- 
formed without any errors or ambiguities 0. However, 
a POVM can be optimal in the sense that it minimizes 
the error probability. Another optimal choice would be 
a POVM that never leads to an error but sometimes to 
inconclusive results, but this choice is not suitable for 
quantum states that are not known exactly. 

The explícit construction of the POVM operator has 
only been possible for a few simple cases. For example, 
if only the two states pA — \A){A\ and ps = \B)(B\ 
can occur with probabiliti es Pa + Pr = 1, an optimal 
POVM can be constructed [15l Í17| via the projectors on 
the eigenvectors of the operator 



A = papa - PbPb 



(9) 



Keeping in mind that the dimension of the Hilbert space 
grows exponentially, finding the corresponding eigenvec- 
tors already becomes a computationally challcnging task. 
More important however is the difnculty that in an ex- 
periment, the measurement of such an observable must 
be implemented. 

Associated with an optimal POVM one obtains the 
minimum error in distinguishing between the states \A) 



6 



and \B) from the Helstrom formula 

Perr = ~ U - ^1 - 4 PA p B |(^|S)| 2 ) . (10) 

Of course, such a simple construction is only possible if 
non-perturbed quantum states are used. Nevertheless, 
one can compare the error probability of the quantum 
algorithm in case of quantum state discrimination with 
the optimal one, see figure [S] Without noise, the error 
probability of the quantum algorithm lies only a few per- 
cent above the lower bound of about 6%. In the presence 
of noise, the quantum algorithm is of course more likely 
to make mistakes, but then it should not be compared 
with the POVM bound since this was derived assuming 
exactly known quantum states. For comparison, we plot- 
ted the performance of the naive POVMs containing the 
projector expectation vàlues |^4)(A| and \B)(B\ (dotted 
lines), respectively, which are clearly inferior to the pre- 
sented quantum algorithm in the presence of noise. 

512x512pixel/0<lkl<40 
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FIG. 5: [Color Online] Error probability of the quantum 
algorithm when applied to state discrimination for equally 
distributed states pa = Pb ~ 1/2. The dashed horizontal 
line corresponds to the minimum error fllOP for unperturbed 
states \A), \B) (as in figure |21 top left panels). The errors 
of the quantum algorithm (large circle symbols) correspond 
to the averaged failure probabilities in figure 0] with the re- 
spective template. For comparison, the error probabilities 
obtained from the expectation vàlues of the projection opera- 
tors |A) (j4| and \B){B\ have been added (small symbols). At 
the price of introducing an inconclusive result (hollow triangle 
symbols) , the probability of obtaining a conclusive but wrong 
result can be lowered even further (filled triangle symbols), 
see section IVIIII Note that the POVM bound does not ap- 
ply in this case and that the decline of this error probability 
with increasing noise is due to the fact that the probability of 
obtaining a conclusive result decreases also. The conditional 
probability of getting the wrong answer provided that the al- 
gorithm arrived at a conclusive result of course increases with 
noise - but is still below the other curves (not shown). 

As will be discussed in the following section, it is pos- 



sible to reduce the error probability even further at the 
expense of introducing an inconclusive result. 



VIII. ALGORITHMIC EXTENSIONS 

Note that the algorithm is vulnerable to translation, 
i.e., it does not recognize an image that has been trans- 
lated in position space compared to the template. To 
cope with translations, it is necessary to use more than 
one reflected photon and to analyzc the two-photon cor- 
relations. After determining the center of mass of the 
image, for example, one could perform the same algo- 
rithm with a template that has been shifted accordingly. 

Since the objective is to shine as few photons as possi- 
ble onto the image, one should try to extract as much in- 
formation as possible from single measurements. Unfor- 
tunately, after a complete measurement of all the qubits 
separately, the full quantum state has been projected 
onto the outcome and no information is left. However, 
in order to accept or reject the hypothesis of a given 
template, we do not need to measure all the qubits. It 
is completely sunicient to ask the question "Are all the 
non-ancilla qubits in the state |0) or not?" Alternatively, 
one could also construct a unitary operator (as with the 
noise-filter operator J\f) that uses another ancilla qubit 
to answer this question. If the answer by measurement 
is "yes", we have accepted the hypothesis (probabilisti- 
cally) and the full quantum state has been projected onto 
the state |0 (na=) ) (g> jo^»') ® |1) <g> |1). However, if the an- 
swer is "no" and we have rejected the hypothesis (again 
probabilistically) , the quantum state has only partially 
been projected onto the high-dimensional sub-space or- 
thogonalto |0( na, ))(S>|Cr n! ") and still contains a significant 
amount of information. 

For example, let us assume that we performed the al- 
gorithm with the template A and that this hypothesis 
has been rejected, i.e., not all of the qubits were in the 
state |0). One possibility for this rejection would be that 
the input state was \B) or some other quantum state 
(rejection correct) or that an input state \A) was falsely 
rejected (either due to strong perturbations or just due 
to the inherent probabilistic nature, see figure EJ. In 
this case, we can partially undo the operations specific 
to template A by applying a Hadamard gate onto each 
qubit again and Grover-rotating back with the template 
A. Now we may switch to the template B and perform 
the inverse Grover rotation with the new template B plus 
the Hadamard gates and measure the outcome. In case 
the overlap between the images A and B is small enough, 
the result of this measurement indicates whether the im- 
age is probably B and the rejection of A was correct or 
whether something clsc happcned (e.g., the rejection was 
erroneous or the image is neither A nor B). 

Thus, if the first template had M t p points and the 

(2) 

second template had M tp points, we could continue the 
quantum algorithm as follows 
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11. apply Hadamard gates on non-ancilla qubits 



12. perform iï = 



l^N x N y /M$ 



with respect to templatc A 



|*i 2 ) =^|*ii) 



Grover iterations 



13. switch to the second template B and perform 



R' 



l^N x N y /M { 



(2) 
tp 



inverse Grover iterations 



with respect to template B 



13/ 



Gb*'|*12> 



14. apply Hadamard gates on non-ancilla qubits 

15. measure final state in computational basis 
|* 15 ) = |0<"-)) ® |0C"«)) ® |i) ® |i) 

The measurement at step 10 of the quantum algorithm 
essentially subtracts the template A from the state. Con- 
sequently, in case of a false rejection, the state amplitudes 
will concentrate on the perturbations of the image and 
the interference patterns resulting from the Fourier noise- 
filtering, see figure |ü] left panels. The resulting state is 
then nearly orthogonal to the second hypothesis in gen- 
eral. In case of a correct rejection of the first hypothesis, 
the state's squared amplitudes will decrease where the 
state overlaps with the first template, see figure right 
panels. Note that here also the different magnitudes of 
the amplitudes due to the different number of points will 
in most cases prohibit a complete removal of these am- 
plitudes. In this case, the overlap between the state and 
the second template will be substantial. 




FIG. 6: [Color Online] Color-coded distribution of the 
squared amplitudes of the state after step 12 of the quantum 
algorithm (rejection of the first hypothesis). Color-coding has 
been chosen as in figurcE]and all images have been perturbed 
by 5% initially. From left to right, the combinations have been 
chosen as follows: (image A L· template A), (image B L· tem- 
plate B), (image A & template B), (image B & template A). 

In this second attempt, the recognition probability for 
a correct hypothesis is naturally a bit smaller as if the sec- 
ond hypothesis had been tested directly, see large versus 
small symbols in figure [7| However, this procedure has 
the advantage, that no further interaction with the image 
is necessary. A further important advantage is that the 
probability of falsely recognizing a template is extremely 
low (large circle symbols in figure 01. This is due to the 



small overlap of the second hypothesis template and the 
falsely rejected state vector (compare again figure [ü] left 
panels) . 

512x512pixel / < Ikl < 40 

bold letters (A : 27416 points, B : 39431 points) 
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X3 0,1 
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■■- 


- ■ state A given 


B rejected 
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■■- 


- ■ state B given 


A rejected 


B tested 


•■- 


- # state A given 
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B tested 


#- 


-# state B given 


B rejected 


A tested 



distortion of given letter 

FIG. 7: [Color Online] Probability for acceptance of the sec- 
ond hypothesis after rejection of the first one. For compar- 
ison, the small hollow symbols (same data as large symbols 
in figure _ corresponding to the acceptance probability if 
the second hypothesis had been tested as the first one - have 
been included. 

With two templates given, one has only three different 
outeomes: The first hypothesis can be accepted (1), the 
first hypothesis can be rejected and the second one is then 
accepted (2), or none of the hypotheses is accepted (3). 
In the last case, given that only \A) and \B) can oceur, we 
know that an error must have happened. Consequently, 
the last outeome must be regarded as inconclusive. With 
the triangle symbols in figure [5] it becomes visible that 
with the extended algorithm, the error probability is sig- 
nificantly lowered, whereas the probability for obtaining 
the inconclusive result rises for strongly perturbed im- 
ages as one would expect. Consequently, in alphabets 
containing only few letters, one can identify the letter 
with a reasonable probability and detect false rejections 
with high probability - with just a single reflected pho- 
ton. Thus, the algorithmic extension is similar to a mix- 
ture of minimum error and inconclusive POVMs - with 
the advantage of being applicable to perturbed quantum 
states. 



IX. DISCUSSION AND SUMMARY 

We have demonstrated that a quantum algorithm 
would be far more effective than a classical computer 
in template recognition - firstly regarding the number of 
incident photons and, secondly, in view of the required 
number of gates. 

In case of a successful run, a single photon will suf- 
fice to identify a pattern (probabilistically) without dis- 
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turbing the sensitive image. In this sense, the quan- 
tum algorithm realizes a destruction-free measurement 
(though probabilistically) , which is similar to the Elitzur- 
Vaidman problem Even when one takes the finite 
success probability of the quantum algorithm into ac- 
count, it still constitutcs a major advantage in contrast 
to classical pattern recognition since a few incident pho- 
tons sufiice. In case of hypothesis rejection, the algorithm 
can be used to test a different hypothesis without neces- 
sitating a further photon interacting with the sensitive 
array. (Interestingly, the probability of falsely recogniz- 
ing an input state as the template is much lower in these 
secondary trials than in the first ones.) 

The number of qubits required to run the quantum 
algorithm on a 2 nx x 2™» image is n x + n y + 2 and the 
number of gates scales as n^,+ny for the Quantum Fourier 
transform, which is also much faster than classical meth- 
ods. The number of necessary Grover iterations will be 

rather small 0{^J 2 n * +n y M^ 1 } = 0{l}, since for reason- 

able images the number of points is comparable to the 
number of pixels. Note that the measurements required 
during the quantum algorithm can also be performed at 
the end. 

However, we have only discussed the algorithm from 
a theoretical point of view, neglecting many obstacles: 
For example, we did not discuss the experimental diffi- 
culties that are to be expected during an experimental 
implemcntation of the required quantum circuit. Apart 
from the quantum computer itself, the realization of 
the quantum-controlled refractors poses serious prob- 
lems. Interestingly, for images of moderate size, the 
proposed quantum algorithm can in principle be realized 
with present-day technology using linear òptics, see the 
Appendix. However, the price one has to pay is an ex- 
ponential scaling of the number of gates etc. Therefore, 
only the first advantage (i.e., only a few incident pho- 
tons are required) of the presented quantum algorithm 
survives in that case. 

Despite of these shortcomings, we hope that the pre- 
sented theoretical discussion of the amazing potential of 
quantum algorithms will further contribute to the theo- 
retical as well as experimental developments in this fas- 
cinating area. 



Note added in proof 

If it is experimentally feasible to arrange to photon 
path is a way such that it interacts with the image several 
times, it is possible to increase the efnciency of the initial 
state preparation (i.e., suppress the probability for the 
photon being absorbed) by exploiting the quantum Zeno 
effect along the lines of P. Kwiat et ai, Phys. Rev. Lett. 
74, 4763, (1995). 
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APPENDIX: LINEAR-OPTICS SETUP 

Instead of a binary representation, the quantum state 
can be encoded directly by the path of a single photon 
in the laboratory. This representation has the advantage 
that the proposed quantum algorithm can be realized 
using linear òptics elements which are (in principle ) al - 
ready available with present-day technology, see also |18| . 
Unfortunately, this advantage goes along with a serious 
drawback: The number of elements grows linearly (plus 
logarithmic corrections) with the number of image pixels, 
i.e., the computational complexity is exponential instead 
of polynomial. Hence, this linear-optics scheme is only 
reasonable for images of moderate size. 




FIG. 8: Preparation of the initial superposition state |s). 
A Hadamard gate on the three-qubit-state |000) is emulated 
by the combination of beamsplitters (crossed boxes) and mir- 
rors (sòlid lines), which distribute the amplitude of the in- 
coming photon uniformly on the one-dimensional array (large 
hatched box, partially absorptive and transmittive) . 
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The initial state \s) = WC™*)^™^) <g> WK)|0(™«)) can 
be generated by means of a series of beam splitters as 
depicted in figure |H1 which distribute the photon ampli- 
tude uniformly over the array and thus replace the quan- 
tum controlled refractors. For a one-dimensional array 
of N = 2™ pixels, this scheme necessitates N — 1 beam 
splitters. 

For this implementation, it is convenient to regard the 
array as being partially absorptive and transmittive. If 
the photon passes the array, the quantum state is au- 
tomatically prepared in a superposition of all points as 
in equation On the other hand, if the photon is 

absorbed, no quantum state is created at all, i.e., the ex- 
istence of a photon corresponds to the measurement of 
the first auxiliary qubit. 

Afterwards, a quantum Fourier transform can be per- 
formed with the setup in figure|Ü| In the binary represen- 
tation k = ^ki2 n ~ l and x = ~^2,Xj2 n ~^ , the quantum 

l 3 
Fourier transform, 

possesses the well-known factorization 

QTT\x) = J2\k)U-^ x 

k=0 £=1 ^ l 

í 

x ]J cxp {2ník e x n+m - f 2- m } . (12) 

m=2 

For example, one observes for n — 3 the following phase 
contributions A$ 





= (" 


^k 3 Xl 


exp {27rifc3 


[x 2 2- 2 + x 3 2 




= (- 




exp {27TÍfc2 


h2- 2 ]} 




= (- 


^k 1 x 3 







By adjusting the beamsplitters such that they automat- 
ically induce phase factors of — 1 if (and only if) tra- 
versed vertically on a straight line, we can generate the 
above phase contributions of ir, i.e., the factors of (— 1). 
The remaining phases (in the above example 37r/4, 7r/2, 
and 7r/4) can be implemented by ordinary phase shiftcrs 
placed in the corresponding paths, cf. figure|ÜJ Note that 
further mirrors could be used to bring the input states 
\x1X2X3) in the order of the computational basis. With- 
out these additional phase shifters, the arrangement in 
figure |ü] would correspond to a Hadamard gate on ev- 
ery qubit H {n) . The setup requires tó"" 1 = O(NlnN) 
beamsplitters, i.e., it does not admit the implementa- 
tion of a quantum Fourier transform in a complexity 
that is only polynomial in the number of qubits but only 
reaches the efficiency of the classical Fast Fourier Trans- 
form (FFT). 

Having calculated the quantum Fourier transform, the 
undesired fc-values can be removed by simply placing 




FIG. 9: [Color Online] Linear-optical setup for a quantum 
Fourier transform on a 8-dimensional Hilbert space corre- 
sponding to an array with 8 pixels. After interaction with 
the array, the quantum state (left-hand side) is described by 
the path of the photon (dashed lines). The quantum Fourier 
transform is performed using beamsplitters (crossed boxes), 
mirrors (red sòlid lines) and phase shifters (green boxes) only. 
The distances between the beam paths are assumed to be a 
múltiple of the photon wavelength, such that different dis- 
tances in perpendicular direction do not contribute to the 
relative phase. The beamsplitters are assumed to act with 
a phase factor of — 1 if traversed vertically on a straight line 
(and without any phase otherwise). The phase shifts (num- 
bers beside boxes) have been adjusted according to equation 
H 1 3p . This scheme only reaches the efficiency of the classical 
Fast Fourier Transform (FFT) . Without the additional phase 
shifts, the scheme implements a Hadamard gate on all qubits. 

photon absorbers into photon paths that correspond to 
the respective fc-values (compare steps 5 and 6 of the 
quantum algorithm). The inverse Fourier transforma- 
tion can be performed in analogy to figure |Ü| The im- 
plementation of the inverse Grover rotations (step 8 of 
the quantum algorithm) is also rather simple by means 
of linear òptics: The oracle, which inverts all states be- 
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longing to the template, can be realized by placing a 
transparent glass template of appropriate thickness (n- 
phase) into the optical path. The remaining inversion 
gate 1— \s)(s\ with respect to the equally weighted super- 
position state | s) can be implemented by a similar phase 
gate sandwiched between two Hadamard gates Tí^ 
The required phase gate has to supplement all states ex- 
cept |0) (which is left unaffected) with a phase factor of 
— 1, whereas the Hadamard gates can be implemented as 
in figure [ü] without the additional phase shifters. 

By placing a single-photon detector in the optical path 
corresponding to the state |0 . . . 0) in the computational 
basis without disturbing the other paths, it is even possi- 
ble to perform the final measuremcnt in such a way that 



the subsequent operations in section IVIIII are possible. 
In summary, the presented quantum algorithm can be 
realized experimentally with a few photons using linear 
òptics elements - but the effort (number of devices) scales 
exponentially instead of polynomially. 

Note that, in case the wavelength of the used pho- 
ton is appropriate, an alternative scheme can be realized 
in analogy to optical filtering: A suitably designed Iens 
translates transversal components of wavenumbers (i.e., 
beam directions) into positions in the focal plane and 
thus acts as a quantum Fourier transform, cf. 20] In this 
case, the complexity is still exponential - but somewhat 
hidden in the Iens construction and the needed accuracy 
ctc. 
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